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The topological structure resulting from the network of contacts between grains (con- 
tact network) is studied for very large samples of disorderly-packed monosized spheres with 
densities ranging from 0.58 to 0.64. The hierarchical organization of such a structure is 
studied by means of a topological map which starts from a given sphere and moves outwards 
in concentric shells through the contact network. We find that the topological density of 
disordered sphere packing is larger than the topological density of equivalent lattice sphere 
packings. 



1 Introduction 

To identify the accessible configurations at the 
grain level and to understand which are the 
possible combinations of such local configu- 
rations which generate the global packing is of 
singular importance, hence it is the necessary 
starting point for any fundamental understand- 
ing of granular matter. Until now the empirical 
investigation of the geometrical structure of 
granular packings have been limited by the 
paucity of accurate experimental data. Indeed, 
after the seminal works of Bernal, Mason and 
Scott (IBernal & Mason I960! IScott 1962ft . it 
has been only very recently that the use of 
tomography has allowed one to 'see' the three 
dimensional structure of such systems and 
explore their geometry from the grain level 
up to the whole packing (ISeidler et al. 20001 
ISederman et al. 20011 IRichard et al. 20031 
lAste et al. 20041 lAste et al. 2005ft . 

In this paper we report and discuss re- 



sults from an empirical investigation by means 
of X-ray Computed Tomography on very 
large samples of disorderly packed mono- 
sized spheres with densities 1 ranging from 
0.58 to 0.64. This study is the largest and 
the most accurate empirical analysis of dis- 
ordered packings at the grain-scale ever per- 
formed. A detailed description of the exper- 
imental methodology and apparatus are re- 
ported in (lAsteet al. 20041 lAste et al. 2005ft . 
We analyzed 6 samples of monosized acrylic 
spheres in cylindrical containers with rough- 
ened walls, and with an inner diameter of 
55 mm, filled to a height of ~ 75 mm. In or- 
der to verify possible effects due to gravity and 
sizes we used two kind of acrylic spheres with 
diameters d= 1.00 mm and and d = 1.59 mm, 
with polydispersities within 0.05 mm. The 



1 The density is defined as the fraction of the volume 
occupied by the balls divided by the total volume of the 
region of the space considered. 
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Figure 1. Average coordination number vs. sample 
density. The filled symbols correspond to the sam- 
ples investigated in the present work. The two sym- 
bols '+' are the values from Bernal 1960; whereas 
the 'x' is from Scott 1962. The line is the best 
lineal - fitting of the A-F data constrained to pass 
through n c = 4 at p = 0.55. 



smaller spheres were used to prepare two 
samples (referred as sample A and C here- 
after) at densities 0.586 and 0.619 containing 
about 103,000 and 143,000 spheres, respec- 
tively. Whereas, the larger spheres were used 
to prepare four samples at densities 0.596, 
0.626 0.630 and 0.64 (samples B, D, E and F) 
each containing about 35,000 beads. A X-ray 
Computed Tomography apparatus (see Sakel- 
lariou et al. iSakellar iou et al. 200 4)) was used 
to image the samples. The two large sam- 
ples (A, C) were analysed by acquiring data 
sets of 2000 3 voxels with a spatial resolu- 
tion 0.03 mm; whereas for the other four 
samples (B, D, E and F) were acquired data 
sets of 1000 3 voxels with a spatial resolution 
0.06 mm. The positions of the center for each 
sphere were retrieved, with a sub-voxel pre- 
cision, by a convolution method applied to 
the se gmented flS heppar d et al. 2004|) dat asets 
(see (lAste et al. 20041 lAste et al. 20051) for 
further details). 



2 Coordination Number 



The average number of spheres in contact 
with any given sphere have been widely 
investigated in the literature of granular 
matter dBernal & Mason 19601 IScott 19621 
ISeidler et al."2~UUUl lAste & Weaire 20001 
ISederman et al. 20011 ISilbert et al. 20021 
IDonev et al. 20041) . Indeed, this is a very sim- 
ple topological quantity which gives important 
information about the local configurations 
and the packing stability and determines the 
cohesion of the material. The average number 
of spheres within a radial distance r from a 
given sphere (n t (r)) increases with the radial 
distance. Its empirical behavior is shown in 
Fig|2| One can note that above r = 0.98c? the 
number of neighbors grows very steeply up to 
a 'knee' at about 1.02c? where a slower growth 
takes place. The definition of the coordination 
number depends therefore on the radial dis- 
tance within which two spheres are considered 
in contact. However, a deconvolution method 
developed in (lAste et al. 20051) can be used 
to distinguish between the contribution from 
the spheres in contact and the contribution 
from other spheres which are near but do 
not touch. We estimate that in the 6 samples 
A-F the average number of spheres in contact 
n c are between 5.81 and 6.97. In fig[T] the 
values of n c are reported as a function of 
the sample densities showing a clear and 
consistent increasing trend with the density. 
Such dependence on the packing density 
(Fig [D nas important conceptual implications 
which are discussed in (Ast e et al. 20 05). 



In this paper we construct 'contact' net- 
works by using different thresholds on the dis- 
tance within which two spheres are consid- 
ered in contact. This produces several different 
contact networks with average coordinations 
(n) that increase with the threshold distance. 
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Figure 2. The average coordination number in- 
creases with the radial distance r within which two 
spheres are considered in contact. 

3 Topological Structure Beyond 
First Neighbors 

Any force-path or any infinitesimal local grain 
displacement must mechanically propagate 
from grain to grain through the network of 
touching grains. The understanding of the 
hierarchical organization of such contact 
network beyond first neighbors is therefore 
crucial. Here we apply to granular matter 
an approach which was originally devel- 
oped for the study of crystalline systems 
dBrunner & Laves 197 It IQ'Keeffe 199Tb! 
Conway & Sloane 1997) and disordered 
foams dAste et al. 19961) . The topological 
structure of crystalline frameworks has been 
intensely studied in terms of the number of 
atoms that are j bonds away from a given 
atom. If we start from a given 'central' atom, 
the first 'shell' (distance j = 1) is made by all 
the atoms in contact with the central one. The 
second shell (distance j = 2) consists of all 
atoms which are neighbors to the atoms in the 
first shell, excluding the central one. Moving 
outward, the atoms at shell j + 1 are all the 
ones which are bonded to atoms in shell j and 
which have not been counted previously. In 
infinite, periodic, crystalline structures with 
no boundaries, the number of atoms per shell 



should increase with the topological distance 
and it has been shown that in several three- 
dimensional crystalline structures the law of 
growth for the number of atoms (Kj) at shell 
j can be described with: Kj = ajj 2 + bjj + Cj 
(with dj, bj and Cj coefficients that might 
vary with j but only within a bounded 
finite interval) (IBrunner & Laves 19711 
IQ'Keeffe 199"Tbl IConway & Sloane 19971 
IGrosse-Kunstleve et al. 1996b . Following the 
definition of O'Keeffe (IQ'Keeffe 199 lab , for 
these crystalline systems, the asymptotic be- 
havior of Kj can be characterized in terms of 
an 'exact topological density' : TD = (ctj) /3 
(IGrosse-Kunstleve et al. 19961) . It has been 
noted that such a topological density is inter- 
estingly related to the geometrical density of 
the corresponding crystalline structure and it 
is a powerful instrument to characterize such 
systems. For instance, it is easy to compute 
that the cubic lattice has Kj = 4j 2 + 2. 
Whereas, spheres packed in a bcc (body 
centered cubic) crystalline arrangement have: 
Kj = 6j 2 + 2 (j > 0). On the other hand 
it has been shown ( |Conway & Sloane"! 997) 
that for Barlow packings of spheres, Kj are 
always in a narrow range within: 

21 7 2 

Wj 2 + 2<K J < (j>0) , (1) 

where the brackets |_- • • J indicate the floor func- 
tion. It has been observed by O'Keeffe and 
Hyde ( |0'Keeffe & Hyde 1996D that for lattice 
sphere packings with coordination number n, 
the general rule holds: Kj — (n — 2)j 2 + 2, im- 
plying therefore a = n — 2. 

Beyond perfect crystalline order very few 
results are known either from theoretical, em- 
pirical or numerical point of view. One can ar- 
gue that Kj must grow with a law comparable 
with the law for a spherical shell: Kj ~ aj 2 ~ 
47rj 2 . However, it is also clear that the shape of 
the growing shell and its roughness can dras- 
tically change the coefficient a (as observed 
in two dimensional cases (lAste et al. 19 96)). 
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Figure 3. Shell occupation numbers vs. topological 
distance. The symbols indicate the different sam- 
ples (as in FiglO and the lines are the best-fits (of 
the growing part only) using the polynomial form: 
Kj = aj 2 + cij + Co- The fits are between j = 2 
and j = 10 (for samples B, D, E, F ) and j = 15 (for 
samples A, C). The data refer to threshold distance 
1.05A 

Moreover, it can be shown (lAste et al. 2 005 ) 
that in some topological networks the law 
of growth can follow an intrinsic dimension 
which is different from the dimension of the 
embedding space (3 in our case). This mecha- 
nism can produce power law growth with ex- 
ponents different from 2, or different behav- 
iors such as exponential - or even faster - laws 
of growth (lAste et al. 2005b . 

We observed that the number of spheres at 
a given topological distance j from a central 
one follows a power law growth (see Fig 13) 
until a critical distance j, above which the 
shells hit the sample boundaries and Kj starts 
to decrease. We verify that a quadratic law 
Kj = aj 2 + Cij + c fits quite accurately the 

observed behaviors of Kj for j < j. This fixes 
the intrinsic dimension for these systems equal 
to 3 (which coincides with the geometrical di- 
mension of the embedding space). 

We find that the coefficient a depends on 
the threshold distance within which neigh- 
boring spheres are considered in contact. In- 
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Figure 4. The coefficient a plotted against the av- 
erage coordination number of the contact network 
(n) shows that disordered packings have larger 
topological densities in comparison with lattice 
sphere packings. 



deed, changes in the threshold distances are 
unavoidably associated with changes in the 
contact network and an enlargement of the 
threshold distance must correspond to a thick- 
ening of the shell. It has been noted previously 
that such a threshold affects also the aver- 
age coordination number n in the contact net- 
work (Fig 01 . In Fig.[4]we show that these two 
quantities are positively correlated: the coef- 
ficient a increases monotonically with n. In- 
terestingly, a comparison with the known be- 
havior in lattice sphere packings (a = n — 2) 
shows that these disordered systems are topo- 
logically more compact than the analogous lat- 
tices. This observation might be relevant when 
the structural stability and rigidity of such 
system are concerned. Understanding whether 
this is a bare consequence of topological disor- 
der or it is associated to the properties of such 
granular systems is a challenging topic which 
must be explored. 

4 Conclusion 

The study of the topological structure of the 
contact network at the grain level shows that 
the average number of spheres in contact in- 
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creases almost linearly with packing density 
and lies between 5.5 and 7.5 in the range of 
densities between 0.58 and 0.64. An extrap- 
olation to the random loose packing density 
(p = 0.55) suggests that at this density the sys- 
tem could have an average number of 4 neigh- 
bors per sphere. This implies the possibility of 
a rigidity percolation transition taking place at 
the random loose packing limit. 

The structure beyond first neighbors, stud- 
ied by means of a topological map, shows that 
the contact network has an intrinsic dimension 
of 3. Surprisingly, we found that the topologi- 
cal density in disordered sphere packings is al- 
ways larger than the topological density in the 
corresponding lattice sphere packings. Such a 
larger topological density is an indication that 
the contact network is more compact in dis- 
ordered systems despite the fact that the geo- 
metrical density is lower. This fact might have 
important implication on the system stabil- 
ity and resilience to perturbations and shocks. 
Notably, this fact is consistent with what was 
observed in dAste et al. 20 05 ) where we found 
that in a certain region around the grains, dis- 
ordered packings can be locally more compact 
than the crystalline ones. 
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